Interference in discrete Wigner functions 
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We analyse some features of the class of discrete Wigner functions that was recently introduced 
by Gibbons et al. to represent quantum states of systems with power-of-prime dimensional Hilbert 
spaces [Phys. Rev. A 70, 062101 (2004)]. We consider "cat" states obtained as coherent superpo- 
sitions of states with positive Wigner function; for such states we show that the oscillations of the 
discrete Wigner function typically spread over the entire discrete phase-space (including the regions 
where the two interfering states are localized). This is a generic property which is in sharp con- 
trast with the usual attributes of Wigner functions that make them useful candidates to display the 
existence of quantum coherence through oscillations. However, it is possible to find subsets of cat 
states with a natural phase-space representation, in which the oscillatory regions remain localized. 
We show that this can be done for interesting families of stabilizer states used in quantum error- 
correcting codes, and illustrate this by analysing the phase-space representation of the five-qubit 
error- correcting code. 

PACS numbers: 03.67.Lx, 03.67.Hk, 03.65.Ca 



I. INTRODUCTION 

Continuous Wigner functions 0, allow for the study 
of quantum systems through the representation of their 
state in phase-space, providing a complete description al- 
ternative to the density matrix. The continuous Wigner 
function W(q,p) is a quasi-probability distribution: it 
is real valued and by integrating it along phase-space 
lines one obtains probability densities of certain observ- 
ablcs (quadratures). However, it cannot be interpreted 
as a true probability distribution because it is not pos- 
itive definite. Negative values are associated with the 
existence of quantum interference, which is itself iden- 
tified as a signal of non-classical behaviour 0, 0] . The 
original definition of the Wigner function has been gen- 
eralized to systems with discrete degrees of freedom 
H H Sll H OH El Among these we find the very rel- 
evant case of quantum computers, which are composed 
by n qubits and therefore have a Hilbert space of dimen- 
sion d — 2™ . The discrete analogues of the Wigner func- 
tion have been recently used to study problems connected 
with quantum information processing such as the phase- 
space representation of quantum algorithms I12J. auan- 
tum state tomography 0, IpI. teleportation |l5l Il6j . 
decoherence in quantum walks |17j . and the phase-space 
representation of quantum error-correcting codes |l8j . 

These discrete versions of the Wigner function have a 
number of peculiar properties which make them rather 
different from their continuous counterparts. On the one 
hand, there is not a unique definition and, according to 
the different proposals, the phase-space can be taken as 
a 2d x 2d grid 12] or a d x d one [TJ]. Here we shall 
point out and study some remarkable features of the dis- 
crete Wigner functions introduced by Wootters and co- 
workers jllL Il3j | for systems with power-of-prime dimen- 
sions, d = p n (we shall focus our analysis on the case 
of qubits, p — 2). These Wigner functions are built on 
a d x d grid with coordinates that belong to the finite 



Galois field GF(d) 0. The use of elements of the finite 
field is the key to the construction of a discrete phase- 
space whose geometrical properties are similar to those of 
the continuous phase-space. This geometrical structure is 
exploited in the definition of W(q,p) and has significant 
consequences that link the Wigner functions to the sta- 
bilizer formalism [2(j and to th e prop erties of complete 
sets of mutually unbiased bases 0,|22,|23|]. As stabilizer 
states and Pauli operators play crucial roles in the very 
definition of this class of Wigner functions, these are ex- 
pected to become a useful tool to analyse some quantum 
information problems, such as those related to stabilizer 
codes, Pauli noise channels, etc |18ll24|. 

Some of our results indicate a potential problem for the 
class of Wigner functions defined by Wootters and co- 
workers. For the sake of clarity we describe this problem 
here in simple terms: these Wigner functions might not 
always be a powerful tool to pictorially display properties 
of quantum states. This is in contrast to the continuous 
Wigner function, which is often helpful to put quantum 
interference on display and to study, for example, its dis- 
appearance due to decoherence 17]. It is known that 
certain states have positive Wigner functions and are of- 
ten referred to (in a loose way) as "classical states" . For 
the discrete Wigner functions such classical states have 
been studied elsewhere 0, HH an d belong to the class of 
stabilizer states. As expected, Wigner functions of states 
which arc coherent superpositions of "classical" states are 
not positive, but display oscillations. In the usual contin- 
uous version of the Wigner function the oscillations are 
localized in phase-space and concentrate in the region 
between the two classical states. The discrete Wigner 
functions we study behave in a drastically different way: 
in general the interference spreads over all phase-space, 
affecting even the regions where the original "classical" 
states are localized. We investigate this property in de- 
tail and find that it is a generic feature. For example, 
we show that as the number of qubits increases the inter- 
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ference tends to spread uniformly over all phase-space. 
Simultaneously, the ratio between the typical value of 
the interference terms and that of the "classical" terms 
decays exponentially. As a consequence, the presence of 
interference becomes harder to identify when using this 
class of Wigner functions. 

However, the discrete Wigner function proposed by 
Wootters and co-workers can be defined in a number of 
ways. Exploiting this freedom we show that given two or- 
thogonal stabilizer (classical) states it is possible to define 
a Wigner function such that all coherent superpositions 
of those states have a natural phase-space representation 
(i.e., one in which the quantum interference is localized). 
This class of states contains some cases that are very in- 
teresting from the quantum information point of view. 
For example, stabilizer quantum error-correcting codes 
[261 l27i |28| define encoded states as arbitrary superposi- 
tions of certain stabilizer states (the logical states). We 
analyse in detail the five-qubit code (that corrects all 
one-qubit errors and encodes one qubit of information 
[2^|) and show how to define a suitable Wigner func- 
tion to represent it. The phase-space representation of 
this class of error-correcting codes is interesting because 
Pauli errors correspond to phase-space translations and, 
therefore, have also a natural interpretation. 

The paper is organized as follows: in Section [n] we 
shortly review the definition of the discrete Wigner func- 
tions. Section IIIII contains a derivation of the expres- 
sions for the Wigner functions of arbitrary superposi- 
tions of two computational states, and a brief analysis 
of some of their features. In Section IIVI we discuss the 
conditions the Wigner function should satisfy to avoid 
the overlap between the interference and the interfering 
computational states. In Section[V]we study the average 
behaviour of the interference fringes. Finally, in Section 
IVII we discuss the phase-space representation of the five- 
qubit error-correcting code. In Section rVIII we summarize 
our results. 



II. DISCRETE WIGNER FUNCTIONS 

We present here a short description of the Wigner 
functions proposed in |ll|. restricting our attention to 
systems with a space of states with dimension d = 2" 
and skipping many details and technicalities that can be 
found in 11, 13. The first step in the construction of 
the discrete Wigner function is the choice of the phase- 
space. This is taken as a d x d grid with "position" and 
"momentum" coordinates (q,p) that are elements of the 
finite field GF(d), which contains a zero element and the 
d — 1 different powers u)i , j — 0, . . . d — 2, of a generating 
element lu. Lines in phase-space are defined as sets of 
points a — (q,p) that solve linear equations of the form 
aq + bp = c, where all elements and operations are in 
GF(d) (for instance, all sums are modulo 2). The field 
structure is responsible for the validity of the following 
geometrical properties of the phase-space: two lines are 



either parallel (i.e., they do not intersect) or they share 
exactly one point; given a line there are d— 1 lines which 
are parallel to it (the d parallel lines form a "striation"); 
every line has exactly d points and there are exactly d+ 1 
different striations. 

The following notational conventions will be useful: 
the symbol A will refer to a line in phase-space; A Q will 
be a line that contains the point a, and \^ a line that 
belongs to the striation k. Accordingly, will denote 
the only line in the striation k that contains the point 
a, and Aq,^ the only line that joins the points a and f3. 
The lines of the form Ao,/j are lines which contain the 
origin, and will be called "rays" ; there is one ray in each 
of the striations. In Figure ^ we show the phase-space 
for a two-qubit system, indicating the five different rays. 
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FIG. 1: The phase-space for a system of two qubits. The 
five different rays are indicated: one is horizontal (p = 0, 
labelled H), one vertical (q — 0, labelled V), and the other 
three oblique (u> k q + p — 0, labelled by the parameter k). 
There is one ray in each striation, and all rays cross the origin. 



The key to the construction of the Wigner function is 
to associate each line A to a rank-one projection operator 
P(A). The Wigner function W{a) is then required to sat- 
isfy certain conditions: first, the sum of its values along 
the line A must be equal to the probability associated to 
the projector P(A): 

J2w(a) = Tr(pP(\)) (1) 

where p is the state of the system. Also, in analogy with 
the continuous case, the inner product of two states must 
be obtained from the Wigner function as: 

Tr(p lP2 ) = dJ2Wi(a)W 2 (a) (2) 

a 

Because of the phase-space geometry, these conditions 
imply that the d lines in a given striation should be as- 
sociated to d projectors corresponding to states in an 
orthonormal basis (which means the Wigner function 
is normalized to one). Moreover, we impose that the 
Wigner function of the "line state" -P(A) is equal to zero 
everywhere except from the line A (where it is constant 
and equal to 1/d). Then, as lines belonging to different 
striations intersect at a single point, the bases associated 
to different striations must be mutually unbiased |2l| . 
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The conditions Q 0) determine that the Wigner func- 
tion at any phase-space point is given by: 



- 1 



(3) 



Thus, W(a) depends on the probabilities for all the states 
related to the lines (one in each striation) which contain 
the point a. As in the continuous case, it is convenient 
to express the Wigner function as the expectation value 
of the so-called "phase-space point-operators" A(a): 



W(a) = Tr (pA(a)) 



A(a) 



(4) 
(5) 



These point-operators are hcrmitian and form an or- 
thonormal basis for the space of operators (using the 
Hilbert-Schmidt inner product). 

To complete the construction above and define the 
Wigner function it is necessary to provide a mapping 
between lines in phase-space and states in Hilbert space. 
We impose over this mapping some additional geometric 
constraints. Let us define a set of d x d unitary operators 
T(q,p) acting on the Hilbert space that faithfully repre- 
sent discrete phase-space translations. By this we mean 
that: i) one such operator is associated to each phase- 
space point, and ii) the product of two such operators 
is, up to a phase, the operator associated to the sum of 
the corresponding phase-space points. It is worth notic- 
ing that this implies that the set of d x d unitaries must 
be closed under the product (up to a phase) and form 
a group. For the association between lines and states 
to respect covariance under translations, the quantum 
state associated to a translated line should be identical 
to the state obtained by acting with the operator T(q,p) 
on the original line state. Also, all the lines in a certain 
striation remain invariant under translations associated 
to the points located in the ray belonging to such stria- 
tion (i.e, the line satisfying the equation aq + bp = c is 
mapped onto itself under the transformation q' = q + qo, 
p' = p+po when aqo + bpo = 0). This implies that the ba- 
sis associated to a certain striation should be formed by 
the common eigenstates of the translation operators as- 
sociated to the points in the ray. For this to be possible, 
of course, such operators must commute. So the phase- 
space construction must be such that the operators in the 
rays define d + 1 mutually unbiased orthogonal bases. 

It is interesting to notice that the idea of representing 
quantum states in the phase-space chosen naturally re- 
quires the splitting of a group oidxd operators into d+ 1 
Abelian subsets defining d + 1 mutually unbiased bases 
in Hilbert space. The elements of the Pauli group are 
natural candidates to be used as translation operators. 
These generalized Pauli operators are defined in terms of 
two binary n-tuples q and p as: 



Pi 



> Qj Pj 



X c 'Z p e t 2 qp . 



(6) 



Here the subscripts indicate on which qubit the opera- 
tors act, the binary n-tuples q, p determine which Pauli 
matrix acts on each of the qubits, and the phase is chosen 
to make the operators hermitian. 

To associate a Pauli operator to each phase-space point 
we must map the n-tuples (q, p) to the field elements 
{q,p). This must be done in a way that ensures the com- 
mutativity of the operators T(q,p) that belong to the 
same ray. The simplest way to find n-tuples satisfying 
this constraint is using matrix representations of the field 
GF(d), i.e. a set of d — 1 non-zero matrices which are 
obtained as the first d — 1 powers of a generating bi- 
nary matrix M (plus the zero matrix). Using such set, 
one can assign the position and momentum n-tuples by 
the following prescription: first, we associate the zero 
n-tuple to the zero element of the field. The n-tuples 
associated to the non-zero field elements can be taken 
as qj — uji^ 1 — > goM J ' — 1 , Pj = uj^ 1 — > PqM^~ x , where 
j = 1, . . . , d— 1, M is the transpose of M and qo, Po are 
arbitrary n-tuples |18| . Given this assignment one can 
show that the Pauli operators located along the same ray 
commute (the Paulis associated to the ray Lo k q + p = 
are generated by the operators T(goM- J ' J pbM J ' + ") for 
j = 1, . . .n). This method is not only a way to build a 
consistent phase-space structure but also a constructive 
way to split the Pauli group into d+1 Abelian subgroups. 
In Figure|21we illustrate the case of two qubits, where M 
is a symmetric matrix and therefore the assignment of 
n-tuples can be chosen to be the same for the position 
and momentum coordinates. 
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FIG. 2: The assignment of n-tuples to the phase-space coor- 
dinates for a system of two qubits; this allows to associate to 
each point a Pauli operator T(q,p) that represents a transla- 
tion in phase-space. 



Once the coordinates are mapped to n-tuples accord- 
ing to this procedure, the association between lines and 
states must be such that the Wigner function exhibits 
covariance under translations: 



W(q,p) 



T(q ,Po) P T^{qo,p ) 
W(q + q ,p + p ) 



(7) 



The method described has an important ambiguity 
(which can be interpreted also as an important degree 
of freedom): it assigns a basis of Hilbert space entirely 
formed by stabilizer states to every phase-space stria- 
tion, but it does not determine which state in the basis 
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we should associate to each line in the striation. As dis- 
cused in [TlT ] , the freedom in this choice leads to a number 
of possible "quantum nets" , each of which is a mapping 
A — * -P(A) between lines in phase-space and states in 
the corresponding bases. The definition of the Wigner 
function will depend on the quantum net chosen, so the 
recipe we have outlined results in a class of different pos- 
sible Wigner functions. At the moment there is no cri- 
terion to choose a unique definition, though there have 
been attempts to reduce the freedom in the choice of the 
quantum net [lH, I3f1 ] . 

Finally, a few notational remarks. It will be conve- 
nient to express the point-operators in terms of the Pauli 
operators in the way |l8j: 



where a A (3 — q a ■ pp — qp ■ p a is the symplectic product 
between a = (q a ,Pa) and (3 — (qp,pp), and the fp are 
constants given by: 

fo = 1, 

/^ = Tr(p(A 0j/3 )T(/3)), (3^0 

Because of the covariance condition, the projectors 
P{\o,p) are eigenstates of the Pauli operators T{j3), and 
so the constants fp are the corresponding eigenvalues. 
Each of them is then equal to ±1, and the choice of the set 
{fp} completely fixes the quantum net. But these con- 
stants are not all independent; for each projector P(Xq^p) 
we can choose only n eigenvalues, corresponding to the 
n Paulis that generate the set of d Paulis on the ray of 
the striation. 



III. INTERFERENCE IN THE DISCRETE 
WIGNER FUNCTION 

Given two states I'i/'i) and \1jj2) with Wigner functions 
W^a) and W^ot), the Wigner function for the superpo- 
sition state \tp) = a\ipi) + b\tp2) is: 

W{a) = \a\ 2 Wi{a) + \b\ 2 W 2 (a) + 

+2TZe{ab*(ij 2 \A(a)\i: 1 )}. (10) 

Thus, the Wigner function of the coherent superposi- 
tion is a weighted sum of the Wigner functions of the 
states \ipi) and IV'2), plus an interference term. If we 
consider states 1^1,2) such that their Wigner function is 
non-negative (line states, for the discrete case discussed 
above |24jp. the interference term is responsible for the 
negativity of the total Wigner function. In the continu- 
ous case, it is well known that the only pure states with 
non-negative Wigner functions are the coherent states, 
whose representations in phase-space are Gaussian. Su- 
perpositions of such states have Wigner functions with 
oscillations which are localized in the phase-space region 



between the two Gaussian peaks. In what follows we 
shall see that for the class of discrete Wigner functions 
with coordinates in GF(d) sketched in the previous sec- 
tion the representation of interference can be drastically 
different. 

We consider a system of n qubits and focus on the 
simplest case: a superposition of two states in the com- 
putational basis. This corresponds to the superposition 
of a pair of vertical lines, and can be easily generalized to 
superpositions of states associated to any pair of parallel 
lines. States in the computational basis are written as: 
\k) = \k\ . . . k n ), with k a binary n-tuple. We shall take 
\il>i) = |0), \ip2) = \'rn) {rn 7^ 0); any other combination of 
computational states can be obtained by translation of 
states of the form considered. The Wigner function W\ 
(W2) takes the value l/rf on the points that belong to 
the vertical line q — (q = m), and vanishes everywhere 
else. The interference term can be calculated from: 

(m\A(q,p)\6) = 1 £ (-l)^'-«"^/(^,p- 

q ',p ' 

(m\T(q\p')\Q)= (11) 

= He*— £W)(-i)* ? eW 

P ' 

According to this expression, the interference term de- 
pends on the position and momentum coordinates in a 
separable form. The only dependence on the momentum 
p is given by an oscillation with "frequency" determined 
by the displacement m between the lines in the super- 
position. The dependence on the position coordinate is 
more involved and is given by: 

F ^ = ^ E W ')(-!)^ ™' P ' (I 2 ) 
v 1 

which depends also on the relative displacement m and 
on the quantum net chosen. 

In contrast with the behaviour of interference in the 
continuous Wigner function, here the interference term 
may be non-zero at any point in phase-space. For any 
fixed pair of states |0) and \fh), it is possible to find 
a quantum net for which the Wigner function mimics 
(as much as possible) the continuous case, being non- 
zero only on a pair of vertical lines. This can be done 
by imposing the condition that the sum l|12|) be to- 
tally constructive at some arbitrary value <f/; then, as 
F(q + m) = F(q)*, it will be constructive at q*i + fh too, 
and it will vanish at any other value of q. 

As an example, let us look at the Bell state \<f>+) = 
— ( 1 00) + \ll))/^/2. This state is a superposition of the 
computational states 1 00) and |11), and is an eigenstate 
of the generalized Pauli operators X^X^) and Z^Z^)- 
This implies that its Wigner function must be invariant 
under the translations in phase-space associated to these 
operators. In [18] it is argued that, taking these symme- 
tries into account, only two possible representations can 
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be obtained for this state; these are shown in Figure |3 
In one of them (Fig a) the Wigner function is nonzero 
only at the corner points; this corresponds to having the 
interference right over the interfering vertical lines. In 
the other representation (Fig|3|b) the Wigner function is 
negative at the central points of phase-space and positive 
everywhere else; this is what we get when the interference 
lies between the original lines. 




00 10 01 11 i 00 10 01 11 



FIG. 3: The possible Wigner functions for the Bell state 
= (|00) + |ll})/\/2, which is a superposition of two com- 
putational states, a) For some choices of the quantum net, the 
Wigner function is nonzero only at the corners of the phase- 
space; this happens because the interference term overlaps 
the interfering lines, b) For the remaining quantum nets, the 
Wigner function is positive at the boundaries of phase-space 
and negative at the central points; this corresponds to having 
the constant original lines and the interference lying between 
them. 

This case shows two qualitatively different ways in 
which interference can behave: it can lie between the lines 
in the superposition, or overlap them; the two possible 
representations look very different from one another. 

IV. THE OVERLAP BETWEEN THE 
INTERFERENCE AND THE ORIGINAL LINES 
CANNOT BE AVOIDED 

In this section we will analyse the conditions the quan- 
tum net should satisfy so that the interference always lies 
outside the interfering vertical lines. This means that 
F(q) i|12|) should be zero at q = (and then it will van- 
ish at q = fh too). This requirement is given by: 

F (0) = ^EW')^ W '-0 (13) 
p ' 

The previous equation can be separated in its real and 
imaginary parts; for each value of fh we then have two 
linear equations on the set f(m,p ')• If we want the inter- 
ference not to overlap the original lines for arbitrary pairs 
of vertical lines, we have a total of 2(d— 1) equations, for 
the d — 1 non-zero values of fh. But the constants f($ t p) 
are not all independent: we can choose arbitrarily the 
values for fig it 0\, where (e*j)j = &y for i,j — 1, ...,n; 
the remaining values are then fixed, and can be obtained 
as products of the independent ones. So the conditions 



on the quantum net are non-linear and become more in- 
volved as the number of qubits is increased. In the follow- 
ing we shall focus on the cases of two and three qubits, 
and see that no quantum net can satisfy the whole set of 
equations for all fh. 

For a system of two qubits, there are three non-zero 
2-tuples fh. For each of them, we can write down the 
conditions ifHfl) as: 

/(io,oo) + /(io,oi) + *(/(io,io) + /(io,ii)) = 
/(oi.oo) +/(oi,io) +*(/(oi,oi) + /(oi,ii)) = (14) 
/(u,oo) - /(ii,u) +*(/(ii,io) + /(n,oi)) = 

The constants f( q nn ~i can be chosen to be 1 without any 
loss of generality pdl|.As mentioned before, the other con- 
stants are not independent: the non-zero points on each 
ray are related to a set of three commuting Paulis gener- 
ated by two of them, so that the third Pauli in each ray 
is (up to a sign) a product of the other two. Taking into 
account the phases of the Paulis T(q,p) © we can see 
that: 

T(ll,10) = -T(10,01) T(01,ll) 
T(ll,01) = -T(10, 11) T(01, 10) (15) 
T(ll,ll) = T(10,10) T(01,01) 

As the f's are the eigenvalues of the translation operators 
©, they must satisfy the same relations: 

/(li.io) = — 7(10,01)7(01,11) 
/(li.oi) = - /(io,ii) /(oi.io) (16) 
7(ii,ii) = /(io,io)/(oi,oi) 

Plugging this in the conditions l|14|) we find that: 
f /(lO.io) + /(oi.oi) = 

\ 7*(10,10)/(01,01) — 1 

which cannot be solved by any choice of /(io,io) an d 
/(oi,oiV Thus, for a system of two qubits there is no 
quantum net that avoids the overlap for arbitrary pairs 
of computational states. 

As the number of qubits in the system is increased, 
the problem becomes much more difficult, not only be- 
cause of the exponential number of equations, but also 
because the maximum degree of the equations is given 
by the number of qubits. For the case of three qubits, we 
wrote the conditions (|13fl and the relations between the 
constants /, obtaining a system of 14 non-linear equa- 
tions for the 21 independent constants. We examined all 
the possible sets of values and none of them solved the 
system. In this way, we found that no quantum net can 
avoid the overlap in the three-qubit case. 

From the results corresponding to systems of two and 
three qubits, it seems a natural conjecture that for higher 
numbers of qubits it will not be possible either to choose 
a quantum net such that the interference lies outside the 
interfering lines for arbitrary superpositions of states in 
the computational basis. 
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V. AVERAGE BEHAVIOUR OF 
INTERFERENCE IN THE WIGNER FUNCTION 

As no quantum net seems to be particularly convenient 
to represent arbitrary superpositions of computational 
states, we studied the "average" representation of inter- 
ference by picking the quantum nets at random. Our goal 
was to analyse its behaviour for large numbers of qubits. 

As was pointed out in section IIIII for each value of rh 
we can always find a quantum net such that the inter- 
ference concentrates in two vertical lines (with relative 
displacement to). In this case, over those two lines the 
real and imaginary parts of F{q) satisfy \Re (F(q)) \ = 
= \Im (F(q)) | = l/2d; this means that the values of the 
interference term are of the same order as those of the 
Wigner function on the interfering lines. However, if we 
choose a different value for to and keep the same quantum 
net, the interference will not retain this nice behaviour, 
and it will spread and fade out. 

We studied the "average" properties of interference 
in the following way: we calculated the distributions 
F(q) for random quantum nets, and we worked sepa- 
rately with the absolute values of the real and imaginary 
parts, defining the new functions R(q) = \lZe (F(q)) |, 
I(q) = \Tm (F(q)) |. We found the maxima for these two 
functions, so as to compare the order of magnitude of 
the interference with the value of the Wigner function 
on the original lines (which is of order 1/d). Besides, 
we calculated the entropy associated to the distributions 
R(q), I(q}- hi order to do so, they were normalized to 
one, so that their values could be interpreted as probabil- 
ities V(q). To facilitate the comparison between systems 
with different dimensions, we defined a normalized en- 
tropy given by: 

S = -Y,V{q) log d {V(q)) (18) 

Q 

Then, the entropy for a distribution concentrated in only 
one value of q is equal to zero, while the maximal entropy 
corresponds to the uniform probability distribution, and 
takes a value of one for any number of qubits. 

Maximal values and entropies for the distributions 
R(q) and I(q) were calculated for 50 random quantum 
nets, and for each possible displacement m, in systems of 
2 to 10 qubits. For each number of qubits, the 50(c? — 1) 
values were averaged. The results for R(q) and I{q) were 
the same, taking into account their mean deviations. We 
found that the ratio between the average of the maximal 
interference values and those of the Wigner function on 
the original lines decays exponentially in the number of 
qubits. For the entropy we found a saturation trend as 
the system is enlarged; this means that for large numbers 
of qubits the interference is spread over the whole phase- 
space taking almost uniform absolute values. As a conse- 
quence of these two results (shown in Figure interfer- 
ence becomes hard to detect using these Wigner functions 
for systems composed by more than a few qubits. 
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FIG. 4: "Average" properties of interference in phase-space. 
On top: decay of the ratio between the maximal interference 
values and those of the Wigner function on the original lines. 
Below: normalized entropy of the interference distribution. 
The results are averaged over 50 random quantum nets and 
the d — 1 non-zero relative displacements m, and the error 
bars are given by the mean deviations. The continuous line 
shows the exponential fit. 



VI. A WIGNER FUNCTION ANALYSIS OF 
THE FIVE-QUBIT ERROR-CORRECTING CODE 

There are important cases in which we are not inter- 
ested in arbitrary superpositions of all possible pairs of 
states in a given basis. For example, in a stabilizer quan- 
tum error-correcting code the encoded states are super- 
positions of a subset of stabilizer states. If a single qubit 
is encoded, a general state is a superposition of a fixed 
pair of stabilizer states. In such case, we can always find 
a suitable definition of the Wigner function, in which all 
the encoded states have simple representations. Below, 
we will show how to implement and use this construction 
to represent the "perfect" five-qubit code, that encodes 
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one qubit of information and corrects all one-qubit errors 




A. The perfect five-qubit code 

This code protects the state of one qubit by encoding 
it in a five-qubit system. The space of "legal" encoded 
states is formed by all the linear combinations of the 
logical states |0)l and (which are the encoded one- 
qubit computational states). The states \0}l and 
are stabilizer states; more precisely, the state |0)l can be 
chosen to be the graph-state p3ll l33 | associated to a five- 
qubit ring. The stabilizer of this graph is generated by 
the Pauli operators Sj that apply X on the j-th vertex 
and Z on all its neighbours, 

S j = z U-i) x U) z U+i) ( 19 ) 

(where indices are cyclic from 1 to 5). The subspace of 
encoded states is formed by the states that are stabilized 
by all the operators of the form SjSk- The graph-state, 
which is stabilized by the operators SjSk and by the op- 
erator S = Si ... S5 — —Xi-u . . . X(5), is chosen as the 
logical state |0)l. The state is then taken to be the 
state which is stabilized by the operators SjSk-, and by 
the operator —S. This means that S plays the role of the 
encoded Pauli operator Z, that distinguishes between the 
two computational states. The encoded computational 
states are connected by the encoded X operator: 

\1) L =Z (1) ...Z (5) \0) L (20) 



qubit). Then we can redefine the translations as: 

n 

T'(q,p) = Yl X T Z T 'e lf93Pj = X^Z^e**** (21) 
3=1 

The Wigner function is now defined as the expectation 
value of the new point operators, given by: 

A » = iE^ 1 )^ T 'W ( 22 ) 



(where the f's are eigenvalues of the redefined transla- 
tions). The construction proceeds just the same way as 
before; each striation in phase-space will be related to a 
stabilizer basis, and the different bases will still be mu- 
tually unbiased, but the basis that corresponds to each 
of the striations might change. To represent the states in 
the perfect five-qubit code it is convenient to choose the 
generators of vertical and horizontal translations as: 

X'j = Sj = Z u _ 1) X u) Z ij+1) (23) 
Z> = Z U) (24) 

Then, the vertical translations are still products of Z 
operators (and so the vertical lines are still the compu- 
tational states) but the horizontal translations are asso- 
ciated to the stabilizers of the graph-state. In this case, 
the horizontal ray can be chosen to be the graph-state 
(the encoded |0) state), and by applying Z operators to 
it we can obtain the states associated to the remaining 
horizontal lines. 



B. A convenient definition for the Wigner function 

To represent this stabilizer code we can use an alter- 
native definition of the Wigner function, that can be 
built following the same procedure as before but redefin- 
ing the translations in Hilbert space. This will lead us 
to a Wigner function with the same geometrical struc- 
ture, but with a particularly simple representation for 
the states we want to study. An essential step in the 
definition of the Wigner function is the identification of 
the translations with the generalized Pauli operators 
In the previous sections we chose to represent the verti- 
cal translations with tensor products of Z Pauli matri- 
ces and the horizontal translations with tensor products 
of X matrices. For this reason, the vertical lines were 
associated to the computational basis while the horizon- 
tal lines corresponded to the eigenstates of all X opera- 
tors. However, we can choose our vertical and horizontal 
translations in a different way, still using Pauli operators. 
In order to maintain the same structure for the discrete 
Wigner function, we only need to find sets of generalized 
Pauli operators Xj and Zj (j = 1, . . . , n) with the same 
commutation-anticommutation relations as the Xtj\ , Ztj\ 
(note that the X'j and i?j may act on more than one 



C. The representation of the encoded states 

As we mentioned above, the Wigner function we have 
just defined is such that the computational state 
has a particularly simple representation: it corresponds 
to the horizontal ray, p = 0. So its Wigner function takes 
the value 1/d on that line, and vanishes everywhere else. 
The other computational state, can be obtained 

from equation (|2U[) . and is therefore associated to the 
horizontal line p — (1, 1, 1, 1, 1) = fh. All the encoded 
states are superpositions of these two "horizontal" states. 

Following steps similar to those in Section ITTTl we can 
write down the Wigner function for the state — 
= a\0) L +b\l) L as: 

K&P) = \ a \ 2w o&P) + \b\ 2 W[(q,p) + W'MP) (25) 
with: 

W^(q,p) = {1/d) 6(0) (26) 
W[(q,p) = (l/d)6(i>-m) (27) 
W'Mp) = (2/d 2 ) (-if™ 

Ke{ab* £Mf%,™) e ~ lf *•*} (28) 
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where the last sum is over all binary n-tuples x. So, the 
Wigner function is a superposition of the horizontal lines 
associated to the states |0) l and plus an interference 
term that oscillates along the horizontal direction and 
depends on the quantum net chosen. 

Now we can choose a quantum net such that the inter- 
ference lies only on two horizontal lines pi and pj + rh, for 
arbitrary pj. All we need to do is impose the conditions: 



f (x,rn) 



x-pi p t^x-m 



(-1) 



x ■ m 
x ■ rh 



(mod 2) 

1 (mod 2) 



(29) 



(note that this can be done because rh is fixed). Then 
the sum over all n-tuples x will be totally constructive at 
p = pi and p = pi + rh, and will vanish for all the other 
values of p. The resulting Wigner function is given by 
lEH , with: 



W'Mp) 



(-1) 



q ■ m. 



{6(£-pi) [Ke(ab*) +lm(ab*)] + 

+ 5{p-pi-m) [Ke{ab*) -lm{ab*)]\{2,Q) 

In Figure [5] we show the representations of the graph- 
state |0)i, and the state |+) L = (|0) L + \1) L )/V2 that 
corresponds to the encoded version of the one-qubit 
eigenstate of X with eigenvalue +1; all the encoded states 
have representations similar to these. 




FIG. 5: The representation of states in the five-qubit code, 
for a convenient definition of the Wigner function, a) The 
encoded computational state \0)l is associated to the hor- 
izontal ray; b) the Wigner function for the state = 
= (|0)l + |1}l)/v / 2 is formed by the lines associated to the 
encoded computational states, plus a pair of oscillating hori- 
zontal lines; all the encoded states will have a similar repre- 
sentation. 



W 6 '(a) are orthogonal when ^ a W a {a)Wl{a) = 0- As 
mentioned before, all the legal states have representa- 
tions that correspond to constant non-negative values at 
the lines p = and p — rh, plus a couple of oscillating 
lines at p = pi and p = pi + rh. The Wigner function 
vanishes for all the other values of p. 

Because of the covariance condition J7J, the occurrence 
of a Pauli error translates the Wigner function in phase- 
space. For the translated state to be orthogonal to all 
encoded states, it suffices to see that the constant lines 
cannot be translated onto the lines p — or p = rh, 
and that the new oscillating lines cannot be placed at 
p = Pi and p = pi + rh. The coincidence of a translated 
constant (oscillating) line at the place of the original os- 
cillating (constant) lines will not affect the orthogonality 
condition, as the sum over an oscillating line vanishes. 
This shows that the orthogonality between the encoded 
and the faulty states is granted if the translations asso- 
ciated to the errors sum to the p coordinate an n-tuple 
different from and rh. Let us now consider how the 
possible one-qubit Pauli errors transform p: 

1. Zij\ error: it is a vertical translation, that maps 
Pj -+ Pj + 1. 

2. error: it is an oblique translation, equivalent 
to SjZ(j_i\Z(j+\)] its vertical component changes 
Pj-l and p j+ x. 

3. Y(j\: it is also an oblique translation, equivalent to 

SjZ(j_i\Z(j-)Z/j + i\, and so it affects Pj-i, Pj and 
Pj+i- 

As we see, Z errors transform only the component of 
p that corresponds to the faulty qubit, while X errors 
change the components associated to the faulty qubit 's 
neighbours, and Y errors transform the components as- 
sociated to the faulty qubit and its neighbours. None of 
these errors leaves p unchanged, and none changes the 
five components of p. Then, the erroneous states are or- 
thogonal to the subspace of encoded states. Besides, by 
comparison of the transformations of p it is easy to see 
that all one-qubit Pauli errors take the state to orthogo- 
nal subspaces, and so the different errors can be identified 
and corrected. Of course, this is no longer true when we 
consider Pauli errors that act on more than one qubit. 
For instance, an error of the form Zu_i)Zu + x) cannot 
be distinguished from an error -XVj), and so this code 
does not allow for the correction of two-qubit errors. 



To conclude this section, let us analyse the represen- 
tation of the state after a one-qubit Pauli error has oc- 
curred. By means of the discrete Wigner function we 
have defined, we can see that all one-qubit Pauli errors 
take the state to subspaces which are orthogonal to one 
another and to the subspace of encoded states. This is 
nothing but the condition for the error to be detectable 
and correctable by the code [2(| [27j ■ According to equa- 
tion J21, two states with Wigner functions W' a (o) and 



VII. CONCLUSIONS 

We studied the representation of interference in phase- 
space using the discrete Wigner functions defined in . 
We showed that, for a coherent superposition of two sta- 
bilizer states, the interference fringes may spread over 
all phase space, even to the regions where the interfer- 
ing states are localized. We analysed the genericity of 
this behaviour and concluded that for systems of two and 
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three qubits there is no way to define the discrete Wigner 
function (by choosing a particular "quantum net") that 
can avoid this overlap. We conjecture that the same re- 
sult will hold for an arbitrary number of qubits. In the 
light of these results, we analysed the average behaviour 
of the interference term in the Wigner function of a co- 
herent superposition of two fixed stabilizer states. We 
did this by picking random quantum nets, and studying 
the ratio between the maximal value of the interference 
term and the value of the Wigner function on the original 
lines; this ratio seems to decay exponentially with the 
number of qubits. Furthermore, the interference tends 
to spread taking absolute values uniformly distributed 
over all phase-space, as the number of qubits is increased. 



These features suggest that the use of Wigner functions 
in this class may not be appropriate to identify the pres- 
ence of interference in phase-space. 

On the other hand, we showed that it is possible to 
define a Wigner function to represent a class of states in 
phase-space that includes the stabilizer quantum error- 
correcting codes. For this purpose it turned out conve- 
nient to define translation operators that are adapted to 
the code. We applied this idea to analyse the phase- 
space representation of the five-qubit error-correcting 
code. Not only we found a simple representation for all 
encoded states, but showed that single qubit errors also 
have a natural interpretation in this scenario, as they are 
phase space translations. 
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